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Abstract 



In this paper, we discuss the possibihties of reahzing the magnetic multi-valued (MMV) recording 
in a magnetic coupled multilayer. The hysteresis loop of a double-layer system is studied analytically, 
and the conditions for achieving the MMV recording are given. The conditions are studied from 
different respects, and the phase diagrams for the anisotropic parameters are given in the end. 
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Many efforts have been devoted to the studies of magnetic multilayers because there have been lots 
of fascinating behaviors displayed in such systems [1-4] . One of the important applications of magnetic 
multilayers in technology is that they can be used as recording media for memory devices. In such 
materials, the hysteresis loop of one domain should be rather rectangular in order that two messages 
can be recorded in the "spin- up" and "spin-down" states, respectively. Recently, much attention has 
been paid to increase the density of the recording media. One of the proposals is to diminish the 
size of the domain. However, the recording density will eventually come to a limit following this 
way, so that one must try to find new approaches. A simple idea is that: if more messages (than 
two) can be recorded in one domain, the recording density will be highly improved even though the 
domain's size remains the same. This is just the idea of the magnetic multi- valued (MMV) recording 
which is believed to be the next strategy of high density recording and has attracted much attention 
from both experimental and theoretical sides. The MMV recording requires that more (than two) 
metastable phases which are stable enough to record messages must exist in the system; therefore, 
the hysteresis loop for such material should contain more (than one) sharp steps. Experimentally, the 
MMV recording was firstly confirmed by the field modulation method on disks of bilayers [3] or island 
on thin layers [4]. However, the theoretical origin is not yet clear. 

More recently, a quantum theory of the coercive force [5] has been established for magnetic systems 
on the basis of some previous works [6-8]. The concept of metastable state was adopted and the 
magnon excitation gap was found to be an appropriate order parameter to monitor the stability of 
the metastable state and to determine the coercive force of the magnetic systems [5]. The quantum 
approach enables one to study the hysteresis behaviors of a magnetic system from a micromagnetic 
view, and some interesting effects in double-film structures had been discussed by this method [5]. 

The present letter is devoted to proposing a theoretical possibility of achieving the MMV record- 
ing in magnetic multilayers with the help of the quantum method. The main idea is to find more 
metastable states in such systems. It is found that the MMV recording can be realized if the interlayer 
coupling and the perpendicular anisotropic constants may satisfy some conditions. 



In this paper, a double-layer system will be investigated analytically. The Hamiltonian can be 
given by: 

^ = E E Im,m'{R,R')Sm{R) • (R') - hj^ ^mW 

m^m' R,R' m,R 

m R 

where the subscripts m,,m! are the number of layers, and R, R' arc the vectors of lattices on the x-y 
plane, /m.m' (R-i R-') are the exchange parameters and only the nearest-neighbor interaction is consid- 
ered. The single-ion anisotropy is the "easy-axis" case {Dj^ > 0), and the "easy-axis" is perpendicular 
to the film. The spins and the anisotropics in different layers are different. It is supposed that Si > S2 
without loosing any generality. 

Following Refs. 5-6, we will introduce the local coordinates (LC) system {xm,ym, Zm}- The 
spin components in the LC system will have the following relations with those in the original one: 
= cos 0m + smOmS^^Sl^ = Sy^{r),S^ = cos 9mS^ - sin 6m S^. In order to study the 
ground state properties and the low-lying spin-wave excitations, one can apply the usual spin-Bose 
transformation such as Holstcin-Primakoff (H-P) [9] or the complete Bose transformations (CBT) [10] 
to the spin operators in the LC system {S^ , , S^}. In a harmonic approximation, the H-P 
transformation and the CBT yield the same results. Then, after the LC transformation and the Bose 
transformation, the Hamiltonian becomes: 

H = Uo + Hi+H2 + ---, (2) 
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H2 can be written in the momentum k space for x-y plane as following: 
H2 = ^ ^F,„,„-(k,0)a+(k)a„,,(k) 

m,m' k 

+ E E^'»."'('^'^)[«"Wa+,(-k) + a,„(k)a„,(-k)], (3) 

m,ni' k 

in which the coefficients Fm,m' and Gm,m' are defined by: 
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+ ^ Sm'Im,m' C0s(6'„ ~ 6'^„) + h COS 6*™ , (4) 
Fm,m' (k, 0) = — -/„i,„i' \/ SmSm' [1 + COs(0„i — 6m' )], (5) 



G,„,m(k,0) = --^2Sm{2S^-l)Dm sin^ (6) 

Gm,™' (k, 6*) = i/m,m' SmS,n' [1 ^ C0s(6'm - 9^' )] • (7) 

Here, 7fc = (1/-^) ea;p(ik • 6) where the summation 6 runs over the Z nearest-neighbors of a given 
site in the x-y plane. 

In a first order approximation, spin configuration {9m} can be obtained by minimizing the ground 
state energy Uq: 6UQ/69m ~ 0, which yield the following equations: 

Im,m'Sm' Sm{9m - 9^) +hsm9m + Dm{2Sm, - I) SVil9m COS 6',„ = TO = 1, 2, • • • . (8) 

Equations above are just the same as the condition of Hi = 0. The harmonic part of Hamiltonian 
can be exactly diagonalized by a generalized Bogolyubov transformation: 

a+(k) = (k)a+(k)+^y,„^„(k)a„(-k) (9) 

n n 

a™(-k) = ^C/„,„(k)a„(-k) + ^y„,„(k)a+(k), (10) 

n n 

SO that wc get finally: 

H^ul, + J2^Mc^l{^hM + --- (11) 

where the magnon excitation energy em(k) in Eq. ( pi] ) and the coefficients (J7m,n,Kn,n) in Eqs. 
(^-(0) can be obtained from the eigenvalues and the eigenvectors of the following matrix 

.^"^^ ^'S^M- (12) 



i 2g{-k) 



The elements of the sub-matrices JF(k) and tJ(k) in matrix 7i(k) are Fm.m'{9,'k) and Gm,m'{9,^) 
defined in Eqs. (^-(0), respectively [5]. 

Following Ref. 5, the miiiinium value of the magnon excitation energy Cm (k) is defined as the gap: 
A(/i) =Mzn[e„(k)]. 

Eqs. (^ may have many solutions corresponding to various possible spin configurations. For every 
solution of Eqs. (^), one can calculate the magnon excitation gap A(ft.) following the method described 
above. According to Ref. 5, if the gap A{h) is positive, the state described by such a solution is a 
metastable one since a variation from this state must cost energy. However, when the gap comes to 
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zero even negative at a field /ic, such a state will no longer be metastable and a transition from this 
state to another metastable one will take place. Thus, in the case that there are many metastable 
states existing in the system, the MMV recording is possible to take place. 

Eq. (H) has two kinds of solutions: the trivial solutions (i.e. 9m = or,7r, to = 1,2) which 
correspond to the aligned spin states; the nontrivial solutions (i.e. 6m ^ or,7r, m ~ 1,2) to the 
canted spin states. Subsequently, we will discuss both the aligned spin states and the canted spin 
states, and discuss which state the system will transit to if the current state is unstable. 

The following notations will be used. /^^^(R, R') — J, Im,m' (Rj R) = / and Dm{'2Sm ^ 1) = Dm- 
The exchange interaction within a layer should be the ferromagnetic type (J > 0). However, both the 
ferromagnetic and the antiferromagnetic types of interlayer exchange coupling will be discussed (i.e. 
/ > or / < ). 

The aligned spin states: 

In such a system, four aligned states are possible. They are illustrated in figure 1. 
For A configuration, we have 6*1 =0,02 = 0. From Eqs. (|)-(0), we obtain 



f IS2+Di+h+JZS2{l~lk) -IVS^ \ ..ON 

^ ' \ -IVSTS^ ISi+D2 + h + JZSi{l-jk) J ^ ' 

and 

a(k) = (14) 

According to Eqs. (|rTl)-(|l2|), we find that the excitation energy em(k) are just the eigenvalues of the 
matrix .F(k). Thus we obtain: 



AA{h) = h+^[Di+b2+ISi+ IS2 - \J {IS2 - ISi + bi- D2)2 + (15) 

From the discussions above, one may find that the system in A configuration will be stable only in 
the case that 

h > hi, (16) 

where 



hi = ]^[-iSi - 1S2 - bi - 1)2 + ^ iiS2 - iSi + bi- b2f + APS1S2] . (17) 

It is similar for B configuration. One may find that the stable region for B configuration is 

hl<h<hl, (18) 



where 



hi = \[IS2 -ISi+b2~bi + ^J {bi +b2- ISi - IS2Y - 4/25i52], (19) 



hi = i[/52 -iSi + b2-bi- ^ {bi + b2- iSi - 1S2Y - 4/251^2]. (20) 

Considering the symmetry between A,B states and C,D states, it is very easy to understand that 
the stable region of C and D states are [— /ij, —h^] and (—00, —h^?^, respectively. 

The canted spin states 

For every trivial solution, non-trivial solutions can be bifurcated from them at some fields. Around 
the bifurcation points, the variations of the angles from the trivial solution should be very small. Thus, 
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it is reasonable to linearize the non-linear equations to study the behaviors of the non-trivial 
solutions around the bifurcation point. Taking A configuration as an example, we have: sin 0m ~ Qm- 
Then around the bifurcation point, equations (||) will be linearized as: 

iS2{ei-e2) + hei + biei = Q (21) 
iSi{e2-ei) + he2 + b2e2 = Q (22) 



which can be rewritten as: 



- 1S2 1S2 \ f e. 



IS^ ~D2-IS^ J \02 J -^{el)- 
The matrix in the left side of the above equation can be diagonalized with the following two eigenvalues: 



= ~ IS2 -bi-b2 + ^{IS2 - ISi +Di- D2Y + 4/251^2], (24) 



A2 = \[-ISi - IS2 -D1-D2- \/{IS2 - I Si + Di- D2Y + APS1S2]. (25) 



Thus, if ft, = Ai or A2, the linearized equations (|21|)-(|22[) or Eq. (|23| ) may have non-zero solutions. 
If we put this solution as a "guessed solution" into the non-linear equations (H), we can finally get a 
non-trivial solution step by step. The non-trivial solutions bifurcated from other trivial solutions can 
be studied similarly. 

Noting Ai > A2, it is easy to understand that the non-trivial solution can not exist when h > Xi. 
It is interesting to find that Ai = ft^. Thus the canted spin state can appear only in the case that the 
aligned spin state is not stable. 

In order to study whether the canted spin state can be stable or not, one should study the magnon 
excitation gap A(/i) in the vicinity of the critical point for the non-trivial solution. Since the 
angles 9i , 62 are very small in this case, it is reasonable to adopt a first order approximation when 
calculating the magnon excitation gap. From Eqs. the elements of the matrix 7i(ftcjk)|k=o 

in a first order approximation can be given as: 

^m.ra' — ^ni.m' ^~ ^-^rri^m' : (2^) 
Grn^m' — ^m,m' ^Gm^m' ^ (2'^) 

where -Fm,m'; ^m,m' havc been defined in Eqs. by setting k = 0. SFm,m' are found to be: 

<5Fi,i = -1/52(01 - 02)' - (ift° + lDi)el (28) 

SF2,2 = -i/5i(0i - 62)' (ift° + lDi)el (29) 

6Fi^2 = SF2,i = -^Vl^SiS2{9i - 02^ (30) 

Since 9i, 6*2 can be substituted by the solution of the linearized equations (pl|)-(p2|) in the critical point 
h'^ in a first order approximation, they must have the following relation: 

— = ^"^^ = (31) 

O2 ^ IS2 + Di + hO^ ^ ' 

Thus all the terms in Eqs. (p8|)-(p0|) can be obtained after extracting a common parameter 6*2 through 
use of Eq. (]3l|). For example, 

SFi,, = [-1/52(1 - - (f + ^DiK^m. (32) 
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So, based on the perturbation theory, the magnon excitation gap A(/i2) of the canted spin state at 
the critical point /ij? in the first order approximation can be presented as: 



-^1,1 + ^2,2 

A(/i°)/^?| must now be a definite value determined by the parameters. The non-trivial spin state 
which is bifurcated from the B configuration can be studied similarly, and the magnon excitation gap 
A(/i^) for such a canted spin states can also be derived following the same procedure. 

The conditions for MMV recording 

We have studied both the aligned spin state and the canted one. In order to realize the MMV 
recording, the four aligned spin states must be overlapping with each other. Thus, it is required that: 

hi < hi < hi (34) 

On the other hand, to be used for recording, the hysteresis loop should be as sharp as possible. 
Otherwise, it may cause difficulty to distinguish two messages. Thus, the canted spin states should 
not exist: 

^{hl) < 0, ^{hl) < 0. (35) 

(^4|)-(^5|) are just the conditions for realizing the MMV recording in a double-layer structure. The 
conditions are the complicated relations between the single- ion anisotropy parameters {Di,D2), the 
exchange interaction parameter (/) and the spins (5*1, iS'2). We will study them from different respects. 

First, we study what is the requirement for the interlayer exchange parameter / if the two magnetic 
layers are determined. The following model will be investigated: 

Model 1: = 3, = 1, D1/D2 = 2.0 
The critical fields /i", hi, hi have been shown together as functions of I / D2 in figure 2. One may find 
that the exchange parameter should satisfy < I < for condition (p^). If the exchange coupling 
is the ferromagnetic case and is very strong (/ > /^), the B and C states can not be stable at all since 
the two magnetic layers are unwilling to antiparallel with each other (fig. 2). In figure 3, A(/i°) is 
shown with respect to I/D2 in order to study the stabilities of the canted spin state bifurcated from 
A configuration. One may find that: only when the interlayer coupling is the antiferromagnetic case 
and is stronger than a critical value (| J| < could this canted spin states appear. If the coupling 

is the ferromagnetic case, this canted spin state is not able to be metastable. A(/i^) has also been 
studied for model 1, and it is always negative. In all, the exchange coupling should not be very strong 
compared to the anisotropy in order to realize the MMV recording. An example has been shown in 
figure 4 where I/D2 = 0.5. The multi-step shape of the hysteresis loop can clearly observed. 

However, there remains a question. At the field hi where A spin configuration is no longer stable, 
B and D spin states are both stable. Why the system will transit to B spin configuration instead of D 
configuration (fig. 4) ? This question can be answered by studying the magnon excitation spectrum. 
According to Eqs. (|9|)-([l2|), one can get the concrete forms of the Bose operators a„i(k). We only 
study the lowest mode of spin wave, so that k = 0. Suppose m = 1 without losing generality, thus 

e^ihl) = A{hl)=0, (36) 
ai = C^i.ioi + [^1,202 (37) 

where 



C/m = (38) 
iIS2 + hl + Di)^+PSiS2 

IS2 + hl + Di 

{iS2 + hi + biY + PS1S2 



U,.2 - ^^- + ^° + ^^ (39) 
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Since the excitation energy of this mode is zero, if there are any kind of fluctuations, the Bosons 
at this mode must be greatly excited without costing energy. The current spin configuration will 
be completely destroyed because of the excitations. Noting a.i is a linear combination of 01,02, the 
quantities |J7i,to'P may be understood as the possibilities of the Bosons in the m'th layer to be excited, 
thus they must can be considered as the possibilities of the spins in the m'th layer to turn flipping. In 
figure 5, the two quantities |t/i,„i'p are shown together as functions of I / D2 for Model 1. One may 
find that in the region where the MMV recording is permitted, |J7i,iP ~ while |C/i,iP ~ 1- Thus, 
at the field where the A configuration is not able to be stable, the spins in 2nd layer are mostly 
likely to turn flipping while those in the 1st layer are not likely to do so. So, in this case, the system 
will transit to B conflguration instead of D conflguration. One may also find that if the interlayer 
exchange is the ferromagnetic case and is very strong (/ » D2), the two quantities will be close. Thus, 
the two magnetic layers are willing to turn flipping together because of the strong interlayer coupling. 
The MMV recording is not able to be realized then. By the way, if there is no coupling between the 
two magnetic layers (7 = 0), we find that |f/i,ip = and |C/i,2p = 1- This is easy to be understood. 
Because the two layers are not coupled, they can be treated independently. In this case, one may find 
that /i° and h"^ are just the coercive forces of the two layers. So, when the external field reaches 
the 2nd layer will turn over while the first layer will not. 



Finally, We study what kinds of materials can be used for MMV recording. Suppose Si = S2 = 
S, D!^ = {2S—l)Dm/\I\S, the phase diagrams for D[ — D'2 plane are given in ferromagnetic case (figure 
6) and in antiferromagnetic case (figure 7), respectively. The two cases arc quite different. To realize 
the MMV recording, the anisotropics for the two materials can not be very close {D[ ~ D^) if the 
exchange is the ferromagnetic case (fig. 6), while there is no such restriction for the antiferromagnetic 
case (fig. 7). However, a common requirement in the two cases is still a weak interlayer exchange 
interaction. 

If the structure is more complicated, one may easily understand that there may be more sharp 

steps in the hysteresis loop. Thus even more messages can be recorded in one domain. The extension 
from the present double-layer system to a complicated one is rather straightforward. 
Conclusion 

In conclusion, we have analytically studied the hysteresis loop of a double-layer magnetic system. 
We find that when the interlayer coupling and the anisotropies of the two materials satisfy some com- 
plicated conditions, more metastable states will be possible to appear, and the magnetic multi-valued 
recording may be realized. The conditions are discussed from different respects, and the permitted 
values of the anisotropies for realizing the MMV recording are presented. 
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Captions: 

Figure 1: Aligned spin configurations in a double- layer magnetic system. 

Figure 2: Critical fields /i^, and as functions of the interlayer coupling constant for the double- 
layer system. 

Figure 3: A(/i°)/0| as the function of the interlayer coupling constant for the double-layer system. 

Figure 4: Hysteresis loop of a double-layer magnetic system with a ferromagnetically interlayer cou- 
pling constant I/D2 = 0.5. 

Figure 5: The possibilities for the two magnetic layers to turn fiipping (|?7i,ip, |f/i,2|^) as functions 
of the interlayer exchange parameter. 

Figure 6: Permitted values of the two anisotropics for realizing the MMV recording in the ferromag- 
netic coupling case. One can achieve the MMV recording in region I, and can not do so in region 

n. 

Figure 7: Permitted values of the two anisotropics for realizing the MMV recording in the antiferro- 
magnetic coupling case. One can achieve the MMV recording in region I, and can not do so in region 

n. 
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